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ABSTRACT 

While galactic bulges may contain no significant dust of their own, the dust within 
galaxy discs can strongly attenuate the light from their embedded bulges. Furthermore, 
such dust inhibits the ability of observationally-determined inclination corrections to 
recover intrinsic (i.e. dust free) galaxy parameters. Using the sophisticated 3D radia- 
tive transfer model of Popescu et al. and Tuffs et al., together with Driver et al.'s recent 
determination of the average face-on opacity in nearby disc galaxies, we provide simple 
equations to correct (observed) disc central surface brightnesses and scalclengths for 
the effects of both inclination and dust in the B, V, /, J and K passband. We then 
collate and homogenise various literature data sets and determine the typical intrinsic 
scalelengths, central surface brightnesses and magnitudes of galaxy discs as a function 
of morphological type. All galaxies have been carefully modelled in their respective 
papers with a Sersic i?^/" bulge plus an exponential disc. Using the bulge magni- 
tude corrections from Driver et al., we additionally derive the average, dust-corrected, 
bulge-to-disc flux ratio as a function of galaxy type. With values typically less than 
1/3, this places somewhat uncomfortable constraints on some current semi-analytic 
simulations. Typical bulge sizes, profile shapes, surface brightnesses and deprojected 
densities are provided. Finally, given the two-component nature of disc galaxies, we 
present luminosity-size and (surface brightness)-sizc diagrams for discs and bulges. We 
also show that the distribution of elliptical galaxies in the luminosity-size diagram is 
not linear but strongly curved. 

Key words: galaxies: fundamental parameters — galaxies: photometry — galaxies: 
spiral — galaxies: structure — ISM: Dust, Extinction — radiative transfer 



1 INTRODUCTION 

Although the bulge-to-total (B/T), or bulge-to-disc {B/D), 
flux ratio was the third, not prime, galaxy- morphology cri- 
teria employed by Sandage (1961), its systematic behaviour 
along the Hubble sequence is well known (e.g. Boroson 1981; 
Kent 1985; Kodaira, Watanabe & Okamura 1986; Simien & 
de Vaucouleurs 1986). While one-third (two-thirds) of the 
stellar mass density in the Universe today is known to re- 
side in bulges (discs) (e.g. Driver et al. 2007), the actual 
allocation of stars as a function of galaxy type is not se- 
curely known. This is however a quantity of interest. If 
(large) bulges, including elliptical galaxies, are the result 
of mergers, while discs formed via the gravitational collapse 
of a rotating proto-galactic cloud, or from the accretion of 
gas around a pre-existing galaxy, then the B/D ratio re- 
flects the dominance of a galaxy's formation mechanisms 



(e.g. Navarro & Benz 1991; Steinmetz & Navarro 2002 and 
references therein). The reliable separation of bulge and disc 
stars is also important because of the connection between su- 
permassive black hole mass and the physical properties of 
the host bulge (see the review in Ferrarese & Ford 2005). 

For many years it was believed that the bulges of disc 
galaxies were universally described by de Vaucouleurs empir- 
ical i?U-i model (de Vaucouleurs 1948, 1958, 1978). While 
deviations were occasionally noted for individual galaxies, 
such as the Milky Way (Kent, Dame & Fazio 1991), the 
above belief only began to change in earnest after An- 
dredakis & Sanders (1994) showed that the exponential 
model provided a better description to the distribution 
of stellar light in the bulges of 34 late-type spiral galax- 
ies; a result reaffirmed by de Jong (1996a) and Courteau, 
de Jong & Broeils (1996) using larger samples. However 
the fuller renaissance, if one can use such a term, started 
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Figure 1. Average logarithm of the B-band, bulge-to-disc flux 
ratio as a function of galaxy morphological T-type. The solid line 
traces data from Simien & de Vaucoulcurs (1986) who used an 
R'^/* bulge + exponential disc model. The i?^/* model is known 
to over-estimate the bulge flux when the bulge is better described 
with a Sersic B,^^" profile having n < 4 (e.g. Trujillo, Graham & 
Caon 2001). The dashed line traces data from Graham (2001a) 
who applied an R^^" bulge + exponential disc model to de Jong 
&; van der Kruit's (1994) roughly face-on {b/a > 0.625) sample of 
galaxies. Both curves, however, still need to be properly adjusted 
for dust attenuation. 

with Andredakis, Peletier & Balcells (1995) who showed 
that Sersic's (1963, 1968) R^^" model was particularly well 
suited to describing the light-profiles of bulges. This re- 
sult was subsequently confirmed by other investigations 
(e.g. lodice, D'Onofrio & Capaccioli 1997, 1999; Seigar & 
James 1998; Khosroshahi, Wadadekar & Kembhavi 2000) 
and Sersic's model has since become the standard approach 
for describing bulges in both lenticular and spiral galaxies 
(e.g. D'Onofrio 2001; Graham 2001a; Mollenhoff & Heidt 
2001; Prieto et al. 2001; Castro-Rodn'guez & Garzon 2003; 
MacArthur et al. 2003; BalceUs et al. 2007; Carollo et al. 
2007). 

Today, bulges are observed to follow a roughly linear 
magnitude-Sersic index relation, with most bulges having 
Sersic indices n < 4 (e.g. Graham 2001a; MacArthur et al. 
2003). Indeed, Balcells et al. (2003, see their figure la) re- 
vealed that even lenticular galaxies typically have values of 
n around 2, a result which has subsequently been echoed 
by others (e.g. Laurikainen et al. 2006). Application of the 
average B/D ratios derived from studies which used R^^* 
bulge -I- exponential disc decomposition introduces a bias 
into the separation of bulge and disc flux when a bulge does 
not have an ii^^* profile (e.g. de Jong 1996b; Trujillo et al. 
2001; Brown et al. 2003). Specifically, given that the bulges 
of most disc galaxies have Sersic indices n < 4, application 
of the R^'''^ model (with its higher central concentration of 
stars and greater tail at large radii) assigns too much of the 
galaxy flux to the bulge. Figure [T] shows the average B-band 
B/D ratios, as a function of galaxy type, from the data in 
Simien & de Vaucouleurs (1986), whose bulge magnitudes 
were obtained using an R^^"^ model. Also shown in Figure [1] 
are the B-band B/D ratios obtained by Graham (2001a) us- 
ing an 7?^/" bulge model. Fi gure[T]is not however the end of 
the story, as neither study addressed the issue of attenuation 
of the bulge flux due to dust in the disc of the galaxies. 

This brings us to our second, and often over-looked, 
concern. While bulges may contain little to no dust of their 
own — and have therefore in the past not had their flux 



corrected for dust — the galaxy discs which effectively cut 
them in half do contain dust. Consequently, particularly at 
optical wavelengths, one actually sees very little of the star 
light from the portion of a bulge on the far side of a disc. 
The extent to which the bulge and disc flux is dimmed de- 
pends on both the observed wavelength and the inclination 
of the (dusty) stellar disc. The reduction to the observed 
magnitude of the bulge and disc components can be as high 
as 2-3 mag in the B-band (e.g. Driver et al. 2007). In spite of 
this, the overwhelming majority of published (optical) bulge 
magnitudes, and thus bulge-to-disc flux measurements, have 
not taken this into account and are therefore in considerable 
error. 

On the other hand, for over the last decade most studies 
have corrected, at least in part, for the influence of dust on 
the discs. This has been accomplished by noting how param- 
eters such as magnitude and surface brightness change with 
viewing angle and adjusting the observed values to those 
which would be observed with a face-on orientation (e.g. 
Valentijn 1990; de Vaucouleurs' et al. 1991; Giovanelli et al. 
1995; TuUy et al. 1998; Graham 2001b; Masters et al. 2003). 
To obtain the intrinsic (dust-free) values of course requires 
an additional step: namely, correcting the face-on galaxy pa- 
rameters for the influence of dust. 

Sophisticated 3D dust /star galaxy models, which in- 
corporate both dumpiness and explicit treatment of vari- 
ous grain compositions and sizes, now exist. The model of 
Popescu et al. (2000) and Tuffs et al. (2004) self-consistently 
explains the UV/optical/FIR/sub-mm emission from galax- 
ies. Indeed, application of their dust model to 10,095 galax- 
ies from the Millennium Galaxy Catalog (e.g. Allen et al. 
2006) perfectly balances the amount of star light absorbed 
by dust in the Universe today with the amount re-radiated 
at infrared and sub- mm wavelengths (Driver et al. 2008). For 
almost a decade these models accounted for the fact that the 
dust is non-uniformly distributed, allowing for high extinc- 
tion in the central regions, intermediate extinction in the 
spiral arms, and semi-transparent interarm regions which 
accounted for the detection of background galaxies. Their 
detailed model not only allows one to accurately correct for 
inclination-dependent extinction, but also for the additional 
extinction which is present when viewing galaxies with face- 
on orientations. 

The one free parameter in their dust model which pro- 
vides the calibration is the central, face-on, B-band opac- 
ity Tg. Fitting a range of simulated galaxies with varying 
integer- values of r^, Mollenhoff, Popescu & Tuffs (2006) 
provided tables to correct observed disc parameters to their 
face-on, dust-free values. Combining their tables with the 
statistically determined average opacity — 3.8 ± 0.7 re- 
ported by Driver et al. (2007), Section [5] of this paper de- 
rives two new equations which can be used to easily correct 
the observed B, V, I, J and K-hand disc scalelengths and 
central surface brightnesses for both inclination and the at- 
tenuating effect of dust, providing intrinsic, face-on, dust- 
free values. The above face-on, central B-band opacity from 
Driver et al. (2007) was obtained by matching the observed 
inclination-attenuation relations for the Millennium Galaxy 
Catalog data with the dust models of Tuffs et al. (2004). 

In Section [3. II we introduce the galaxy data sets which 
shall be used in our analysis of disc galaxy structural param- 
eters, while Section [3. 2l describes the methodologies adopted 
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to bring this data onto a uniform system. Section |4] provides 
tables of the mean intrinsic structural parameters as a func- 
tion of galaxy type in various passbands. This encompasses 
disc scalelengths, central surface brightnesses and magni- 
tudes. Average K-h&nd bulge effective radii, effective surface 
brightnesses, Sersic indices and magnitudes are also listed. 
In addition, and also as a function of galaxy type, we provide 
the median bulge-to-disc size ratio in the /('-band and the 
median dust-corrected bulge-to-disc flux ratio in the near- 
infrared and various optical bands. 

Finally, a selection of bivariate plots are shown in Sec- 
tion [S] In addition to (surface brightness)-size diagrams for 
discs and bulges, we provide new expressions for the size- 
luminosity relations of discs in early- to mid-type disc galax- 
ies. We also present the size-luminosity relation for bulges 
and elliptical galaxies, emphasizing that it is not a linear 
relation. 

While inclination corrections have been applied for 
many years, and are still necessary today (e.g. Bailin & 
Harris 2008; Mailer et al. 2008; Unterborn & Ryden 2008), 
they only remove one of several systematic biases that cause 
the observed flux distribution to misrepresent the true (in- 
trinsic) stellar distribution of galaxies. Allowances for non- 
homology through the use of Sersic (bulge) models and sep- 
arate dust corrections to both the disc and bulge are vital if 
we are to know the intrinsic physical properties of galaxies. 
It is hoped that the distributions and trends shown here, 
which have been acquired after dealing with the above three 
issues, as done in Driver et al. (2007), will provide valuable 
constraints for simulations of galaxy formation which are 
used to aid our understanding of galaxy evolution (e.g. Cole 
et al. 2006; Springel & Hernquist 2003; Almeida et al. 2008; 
Croft et al. 2008, and references therein). 

A Hubble constant of 73 km s"^ Mpc~^ has been used 
when _ffo-independent distances were not available. 

2 CORRECTIONS FOR DUST AND 
INCLINATION 

2.1 Disc scalelengths and central surface 
brightnesses 

The discs of disc galaxies contain dust, in particular various 
mixes of graphite and silicate. This reduces the amount of 
(ultraviolet, optical and near-infrared) light which escapes 
from such galaxies. The distribution of this dust is known 
not to be uniform, with the centres of discs containing more 
light-absorbing particles than their outskirts (Boissier et al. 
2004; Popescu et al. 2005). Such dust not only decreases the 
observed central surface brightnesses, ^o, of the discs, but 
the radial gradients in the opacity also make the observed 
disc scalelengths, h, greater than their intrinsic values. 

The sophisticated dust/star galaxy model of Popescu 
et al. (2000) and Tuffs et al. (2004) incorporates dumpiness 
and an explicit treatment of various grain compositions and 
sizes within a 3D distribution. Solving the radiative transfer 
equations they are able to measure the radial dependence of 
attenuation in various wavebands. MoUenhoff et al. (2006) 
has used this to provide figures and tables for the effects of 
varying opacity and inclination on the parameters /io and h. 
This was done using a range of galaxy models with different 
integer values for the central, face-on B-band opacity Tg. 



Table 1. Dust correction parameters for equations[T]and[2] which 
have themselves been fitted to the simulated data in Figure |2] 
These equations can be used to correct the observed scalelengths 
and central surface brightnesses of discs for the eff'ects of dust and 
varying inclination, reproducing the stellar distributions' intrinsic 
(i.e. dust-free) face-on values. 
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0.06 ±0.01 


0.79 ±0.02 


1.02 ± 0.01 


0.13 ±0.02 



From a study of 10,095 nearby (0.013 < z < 0.18) 
galajcies in the Millennium Galaxy Catalogue, Driver et al. 
(2007) have recently constrained this mear0 opacity to be 
Tg — 3.8 ± 0.7. We explicitly note that this is a statistical 
determination of the opacity and there may well be galaxy- 
to-galaxy variations outside of these bounds. It should also 
be noted that this value does not make spiral galaxies opti- 
cally thick in their outer parts or interarm regions (Popescu 
& Tuffs 2007). Such high central opacities are however re- 
quired, or rather consistent with, the observed amount of 
infrared and sub-millimetre flux coming from the thermally 
heated dust in galaxies (Driver et al. 2008). Using the above 
value for the opacity, and a linear extrapolation between the 
gridded data points in MoUenhoff et al., the points in Fig- 
ure [5] show the combined effect of dust and viewing angle 
(i.e. inclination) on the disc scalelength and central surface 
brightnesfl in the B,V, I, J and K passbands. The shaded 
region shows the range in behaviour at fixed inclination as 
Tg changes from 3.1 to 4.5. 

Using the linear regression routine FITEXY from Press 
et al. (1992)0 two empirical relations have been fit to the 
simulated data points in Figure [2] They are such that 

MO.oba - AiO,intrin = Oa + [2.5 log(cOS i)] , (1) 

and 

/lobs/ftintrin = f'A - dx log(cOS i) , (2) 

for some inclination i, with i — degrees describing a 
face-on orientation. To assist with the understanding of 
Figure (2] we note that in the absence of dust, the line- 
of-sight depth through a disc increases as the disc is in- 
clined, such that the pathlength changes by (cosi)^^ to give 
Ato,obs = MOantrin + 2.5 log(cos j). The parameters ax, bx, cx 
and dx in equations [T] and [5] are a function of wavelength, 
and given in Table [T] These equations can be used to correct 
observed values to the intrinsic face-on (dust-free) values. 

^ Possible variations of opacity with Hubble Type, and thus also 
luminosity, are not yet firmly established. While Driver et al. 
(2007) and Shao et al. (2007) did not find evidence for such a 
variation, this is still an open issue (e.g. Masters et al. 2003) and 
needs to be addressed when far-infrared and sub-mm data become 
available for galaxies selected from optical surveys. 
^ We have used columns 4 and 5 from the tables in MoUenhoff et 
al. (2006). 

^ The uncertainty on the independent variable, namely the incli- 
nation, was set to zero. The uncertainty on the dependent variable 
(Mo.obs or /lobs) arose from the uncertainty of ±0.7 on the value 
= 3.8 (Driver et al. 2007). 
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Figure 2. Inclination-attenuation corrections to the scalelcngth (h) and central surface brightness (fio) of simulated disc galaxies. The 
data points and the associated shaded area have been extracted from MoUenhoff et al. (2006) using the statistically determined, face-on 
_B-band opacity Tg = 3.8 it 0.7 (Driver et al. 2007). The data points show the difference between the observed values at some inclination, 
i, and the intrinsic face-on (i = 0) values which would be observed in the absence of dust. The solid line is an empirical fit (see equations[T] 
and© using the code FITEXY from Press et al. (1992). The parameters of the fit are given in Table [T1 From 0.7 < 1 — cos^ < 1.0, 
the curve is an extrapolation of the model that was fitted to the data at smaller inclinations; as such it may not be applicable over this 
high-inclination range. 



However, we caution that because the data in MollenhofT et 
al. only extends to inclinations where 1 — cosi < 0.7, the 
curves shown in Figure [2] at higher inclinations are an ex- 
trapolation of these empirical models and may therefore not 
be reliable for systems more inclined than 73 degrees. 

2.2 Bulge and disc magnitudes 

In addition to the above expressions, we shall also use the 
equations in Driver et al. (2008) to correct the observed disc 
and bulge magnitudes for the effects of inclination and dust. 
These are given by 

Mdisc.obs - Mdisc.intrin = bl + &2 [1 - COs(i)]''''' , (3) 

and 

Mbulgcobs - Mbulgo.intrin = dl + d2 [1 - COs(i)]'*^ . (4) 

where the coefficients in the above equations depend on 
wavelength and are provided in Driver et al. (2008) for vari- 
ous passbands, including the B, r' , i' and K bands. We have 
used their r' and i' coefficients for our R and I band data. 
This approximation will not introduce any noticeable bias, 
as inspection of Figure 2 from Driver et al. (2008) will reveal. 

The effect of dust on bulges is far more severe than 
generally realised. Past studies of the B/T ratio, while usu- 
ally correcting the observed disc flux for inclination and 
dust (but typically overlooking any attenuation once ad- 
justed to a face-on orientation), normally assume that the 
bulge is not affected by dust. This is wrong. Moreover, 



given that bulges and discs have their own inclination- 
attenuation correction, which is not surprising, there is no 
single galaxy inclination-attenuation correction. That is, one 
must perform a bulge/disc separation. Unfortunately this 
may severely hamper the ability of large-scale studies (which 
have fltted single Sersic models to two-component disc galax- 
ies) to recover intrinsic galaxy parameters. 

In this study we use the above formula to provide dust- 
corrected bulge magnitudes and dust-corrected bugle-to-disc 
flux ratios. However the influence of dust on the individ- 
ual Sersic parameters of the bulge {fj.c,Rc and n, see equa- 
tionO is not well known. It will of course depend not only on 
the wavelength of one's observations but also on the precise 
star/dust geometry, with concentrated high-n bulges likely 
experiencing a greater degree of attenuation. Moreover, the 
overwhelming majority of optical bulge parameters reported 
in the literature are dependent, to some degree, on the hap- 
penstance inclination of their galaxy's disc. As yet unknown 
corrections to a dust-free configuration are required if we 
are to have an accurate set of physical parameters for the 
bulges. Nonetheless, for now, we have tabulated K-hand 
Sersic bulge parameters as these will be the least affected 
by dust. 



3 THE GALAXY DATA 

In this study we primarily investigate the properties of 'late- 
type galaxies', and more generically 'disc galaxies'. Today, 
the term 'early-type galaxy' is often used to refer to both 
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elliptical galaxies and lenticular galaxies. While the average 
properties of elliptical galaxies are not considered to be a 
function of inclination — and so one should identify and ex- 
clude such systems from current inclination-dependent anal- 
yses of the effect of dust — this is not true for systems 
with large scale discs. We have therefore included lenticular 
disc galaxies in our analysis, rather than grouping/excluding 
them with the elliptical galaxies in an early-type galaxy bin. 

However we note that elliptical galaxies are certainly 
not always devoid of dust (e.g. Ebneter & Balick 1985; Sadler 
& Gerhard 1985; Ebneter, Djorgovski & Davis 1988; Leeuw 
et al. 2008). Indeed, they also frequently possess small nu- 
clear dust lanes and dusty nuclear discs (e.g. Ferrari et al. 
1999; Rest et al. 2001) which probably originate from sources 
such as type la SNe and stellar mass loss from, for example, 
the winds of red giant stars (Ciotti et al. 1991; Calura et al. 
2008). On a galactic scale, however, the dust is destroyed 
rather quickly by ion sputtering from the hot X-ray gas 
(e.g. Temi, Brighenti & Mathews 2007). In contrast, lentic- 
ular galaxies, such as the Sombrero galaxy, are commonly 
observed to possess considerable amounts of dust in their 
discs. In this regard lenticular galaxies are more like spiral 
galaxies. Interestingly, Bekki, Couch & Yasuhiro (2002) and 
Aragon-Salamanca, Bedregal & Merrifield (2006) have ar- 
gued that lenticular galaxies may be spiral galaxies in which 
the star formation has been turned off and the spiral pat- 
tern dispersed. If such a morphological transformation oc- 
curs, one would simultaneously require the bulge-to-disc flux 
ratio to increase. This may proceed via the passive fading of 
the disc and, as noted by Driver et al. (2008), upon removal 
of centrally located dust. 



3.1 Literature catalogues 

There has been a number of papers which have provided 
structural parameters for nearby disc galaxies. At the end 
of 2005, when selecting which of these would be suitable 
for our compilation of bulge and disc parameters, we ap- 
plied the following criteria: (i) a Sersic bulge model had 
been simultaneously fit with an exponential disc model; (ii) 
roughly a couple of dozen or more galaxies had been mod- 
elled; (iii) morphological (Hubble or T)3 types existed for 
the galaxies; (iv) the galaxies appeared large enough that 
their bulges were well-resolved (this meant that we only used 
studies with galaxy redshifts typically less than 0.03-0.04); 
(v) studies of early-type galaxies (i.e. E and SO) which had 
not checked if a bulge-only fit was superior to a bulge plus 
disc fit were excluded as many of the fitted discs (and hence 
the structural parameters) may be spurious; and (vi) stud- 
ies in which the fitted disc frequently failed to coincide with 
the outer light-profile were also excluded, which included 
studies which often had the fitted bulge model contributing 
more fiux than the disc model at the outer light profile. This 
latter problem is illustrated in Graham (2001a, his figure 7) 
and results in overly large Sersic indices, B/T flux ratios 



* We used de Vaucoulcurs' T-typc classification from the RC3 
catalogue (de Vaucoulcurs et al. 1991) which is roughly such that: 
-3 < T < -1 (SO), T=0 (SO/a), T=l (Sa), T=2 (Sab) ... T=9 
(Sm), T=10 (Irr). 



Table 2. Number of galaxies in each passband from the Hter- 
ature data that we were able to use. G02 = Graham (2002), 
low surface brightness (LSB) galaxies only, supplemented with 
LSB galaxies from Graham & do Blok (2001) for which redshifts 
have since become available; M04 = MoUcnhoff (2004); MCH03 
= MacArthur, Courtcau & Holtzman (2003); GOl = Graham 
(2001a, 2003); MHOl = MoUenhoff & Heidt (2001); SJ98 = 
Scigar & James (1998); BGP03 = Balcells, Graham & Pclctier 
(2003); KdJSBOa = Knapcn ct al. (2003); GPP04 = Grosb0l, 
Patsis & Pompei (2004); LSB05 = Laurikainen, Salo & Buta 
(2005); DD06 = Dong & De Robertis (2006); KdJP06 = Kassin, 
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and Rc/h size ratios. This can happen when using a (signal- 
to-noise)-weighted bulge + disc fitting routine on a galaxy 
which has additional nuclear components and, as noted by 
Laurikainen et al. (2007) and Weinzirl et al. (2008), when 
significant bars are present. Our final sample selection, while 
almost certainly not all inclusive due to unintentionally over- 
looked references, is given in Table (2] Unfortunately the op- 
tical data from Hernandez- Toledo, Zendejas-Dominguez & 
Avila-Reese (2007) and Reese et al. (2007) appeared too 
late for us to include. 

Individual studies meeting the above criteria contained 
on the order of ~20 to ~100 galaxies. Previously, when these 
galaxies had been binned into their morphological type (e.g. 
Sab, Scd, etc.), this tended to result in less than ideal num- 
bers per bin in each of these individual studies. As a result, 
past answers to questions of a statistical nature, such as the 
typical value of, and range in, some physical parameter for 
a given morphological type, were poorly defined. By com- 
bining here (see Section |3]2]) the data from these studies, we 
hope to better answer such questions. Importantly, we also 
apply corrections to the magnitudes of both the discs and 
the bulges (see Section [Sjl. 

We stress that while the galaxies used here are typi- 
cal representations of galaxies along the Hubble sequence, 
we have applied no sample selection criteria to provide a 
magnitude-limited sample. Given the inclination-dependent 
reduction to magnitudes due to dust, we note in passing that 
such a task is not as simple as one may first think: after cor- 
recting for dust, one's original magnitude- limited sample will 
subsequently have a rather jagged boundary. However the 
desire to include as many galaxies as possible, rather than 
this issue, has been the driving ratioanle employed here. 
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3.2 Homogenisation 

We have started with the catalogues of Sersic-bulge and 
exponential-disc parameters from the studies listed in Ta- 
ble [2] While the exponential model used to describe the ra- 
dial stellar distribution in discs has been around for a long 
time (e.g. Patterson 1940; de Vaucouleurs 1957; Freeman 
1970), Sersic's (1963) 3-parameter model has only become 
fashionable over the last decade. It can be written as 

fiiR) = Mc + 1.0866„ [{R/Rcf/" - l] , (5) 

where ~ —2.5 logic is the surface brightness at the effec- 
tive half-light radius Rc, and n is the Sersic index quantify- 
ing the radial concentration of the stellar distribution. The 
quantity fo„ is not a parameter but instead a function of n 
such that r(2n) = 27(271,6^), where F and 7 are the com- 
plete and incomplete gamma functions as given in Graham 
& Driver (2005). For 0.5 < n < 10, 6„ « 1.9992n - 0.3271 
(Capaccioli 1989). The exponential model used to describe 
discs can be obtained by setting n = 1 in Sersic's model, 
however it is more commonly expressed as 



fi{R) = fio + im6{R/h), 



(6) 



in which fj,o is the central {R — 0) surface brightness and h 
is the exponential disc scalelength. 

The apparent magnitude of an exponential disc is such 

that 

mdisc = Mo - 2.51og(27r/i^), (7) 
and that of a Sersic bulge is given by the expression 



mbuigc = - 2.5 log(27ri?e) - 2.5 log 



(8) 



with both h and Re in arcseconds, and fio and Hc in 
mag arcsec"'^. The gamma function F has the property 
2nF(2n) — r{2n + 1) = (2n)!. Conversion from apparent 
magnitudes m to absoulte magnitudes M proceeds via the 
standard expression m — M = 25-|-51og(Distance[Mpc]). Al- 
though not done here, applying one's preferred stellar mass- 
to-light ratio to the absolute luminosities yields the associ- 
ated stellar masses. 

In practice, because we use the observed disc magni- 
tudes, due to the apparent ellipticity e of these discs — as 
they are seen in projection - equation [7] had the term 
replaced with h^{l — e). The quantity 1 — e is equal to the 
observed minor-to-major axis ratio of the disc. As we do not 
have information on the ellipticity /triaxiality of the bulges, 
this work has assumed they are spherical (but see Mendez- 
Abreu et al. 2008). Unless specified in their respective pa- 
pers, the inclination of each galaxy has been estimated using 
cos i — b/a, where b/a is the observed outer axis ratio. Discs 
of course have a finite thickness, and so the applicability of 
this expression deteriorates when dealing with increasingly 
edge-on galaxies. In the past, discs have been modelled as 
transparent oblate spheroids (e.g. Holmberg 1946; Haynes 
& Giovanelli 1984; Guthrie 1992) giving rise to the relation 



= [{b/o 



Q ]/[l — Q ], where Q = c/a is the in- 



trinsic short-to-long axis ratio. Transparent edge on discs 
will have an ellipticity equal to Q rather than zero. Not sur- 
prisingly, the above relation breaks down in the presence of 
obscuring dust (MoUenhoff et al. 2006. their section 4.3 and 
Fig. 12 which plots [b/a]obs/[cos j]truo). 



The recovery of disc inclinations is not a major problem 
here because the majority of our galaxies have observed Eixis 
ratios b/a > 0.34 and therefore inclinations less than ~70 
degrees, where the use of cos i — b/a remains a reasonable 
approximation. Only one galaxy (UGC 728, b/a = 0.284) 
in our optical data has b/a smaller than 0.34. In our K- 
band sample of over 400 galaxies, where the effect of dust 
on magnitudes and scalelengths are small, there are three 
samples with (some) discs having b/a < 0.34. These include 

(1) Dong & De Robertis (2006) where a mere six per cent 
of the galaxies have inclinations greater than 75 degrees and 

(2) MoUenhoff & Heidt (2001) from which only six galaxies 
have b/a < 0.34. As such, these few galaxies have no sta- 
tistically significant impact on the results. However, for the 
19 early-type disc galaxies from (3) Balcells et al. (2003), 
the observed ellipticities reach as high as 0.83. Unfortu- 
nately MoUenhoff et al. (2006) do not provide corrections 
to cos i for inclinations greater than 73 degrees, and so we 
have adopted the /C-band value for Q — c/a = 0.11 which 
was used in the dust model of Tuffs et al. (2004)0 While 
this is not ideal, although dust is less of an issue in the K- 
band, we note that for the most edge-on disc galaxy, with 
an observed ellipticity equal to 0.83, the derived inclination 
changes from ~80 degrees to 82.5 degrees, i.e. cos i only 
changes from 0.17 to 0.13. 

From the papers listed in Table [2] we took, when 
available or derivable, the five main structural parameters 
{fj.0, h, fic,Rc,n). We also took the inclination i of the disc, 
which, as noted above, usually came from the ellipticity e 
of the outer isophotes such that e = l — &/a = l — cosi). 
Potential intrinsic disc ellipticities, i.e. non-circular shapes 
of face-on discs (Rix & Zaritsky 1995; Andersen et al. 2001; 
Barnes & Sellwood 2003; Padilla & Strauss 2008), and lop- 
sidedness (e.g. Kornreich et al. 1998; Bournaud et al. 2006; 
Reichard et al. 2008), have not been measured for our galaxy 
sample and are consequently ignored. However given the in- 
trinsic disc ellipticity has been estimated to have a mean 
value of only ~0.05, the uncertainty this introduces to the 
corrected scalelengths via equation [5] is less than 1 percent, 
and the uncertainty on the central surface brightnesses ob- 
tained via equation[T]is less than 0.04 mag arcsec^^. Except 
in the K-hand, one can see from Figure [5] that the uncer- 
tainty in the face-on opacity will introduces a greater source 
of scatter. 

For the distances to the galaxies we consulted Tonry 
et al. (2001) in the case of the lenticular galaxies, and used 
the (Virgo + GA + Shapley) distances given in NEE0, for 
the remainders. As a consequence, a Hubble constant of 73 
km s~^ Mpc~^ has effectively been assumed. This Hubble 
constant is also the value reported in (Blakeslee et al. 2002) 
and is the halfway point between the two values reported in 
van Leeuwen et al. (2007). 



^ The average _B-band value of Q from Xilouris et al. (1999), 
which was used by Tuffs et al. (2004), is 0.074. The intrinsic 
stellar scaleheight reported by Xilouris et al. (1999) is largely 
independent of wavelength, while the intrinsic stellar scalelength 
decreases with wavelength. The dust model of Tuffs et al. (2004) 
used a K-hand intrinsic stellar scalelength that was 0.683 times 
the -B-band value, giving rise to the value 0.074/0.683 ~ 0.11. 
^ NASA Extragalactic Database (NED, 

|http: / /nedwww.ipac.caltech.edu| . 
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Inclination-corrections to the disc surface brightnesses 
which had been applied in some papers were undone accord- 
ing the prescription in each paper. The resultant 'observed' 
disc surface brightnesses were then corrected to a face-on, 
dust- free value using equation [1] Scalelengths — which had 
not been adjusted for inclination — were corrected using 
equation [2l and disc magnitudes were subsequently com- 
puted using equation [7] 

Disc magnitudes were also derived using a second ap- 
proach, in which the observed (uncorrected) magnitude was 
corrected using equation [S] Similarly, the observed bulge 
magnitudes have been corrected here using equation Pub- 
lished Sersic bulge parameters are not corrected for dust; ex- 
pressions to do so do not exist. We therefore report on only 
the if -band Sersic bulge parameter^ which are summarised 
in Table 

Galactic extinction corrections from Schlegel et al. 
(1998), as listed in NED, have been applied (if not already 
done so in the original paper). Cosmological redshift dim- 
ming also reduces the galaxy flux slightly, so we have applied 
an adjustment of 2.5 log(l -I- z)* to the magnitudes and sur- 
face brightnesses when not already done so. This is not a 
huge adjustment, for example, it amounts to 0.1 mag at a 
redshift corresponding to 7,000 km s~^. Given the proximity 
of the galaxy samples, we have not applied evolutionary nor 
Jf-corrections. 

In what follows we have grouped all SO galaxies (—3 < 
T < — 1) into a single bin, denoted by T = —1 in our Figures 
and Tables. 



4 STRUCTURAL PARAMETERS AND THEIR 
RATIOS 

4.1 The near-infrared 

We have computed the median (and width of the central 68 
per cent) from the distributions of various structural param- 
eters as a function of galaxy type. In Table [3] we report the 
face-on, dust-free, _R'-band disc central surface brightness no 
and scalelength h, obtained using equations [T] and [2] These 
values of no and h have been used to derive, via equation [3 
the disc magnitudes given in Table [S] they are also shown in 
Figure [3] Spiral galaxies are commonly referred to as 'early' 
type or 'late' type (Hubble 1926)|!|. Given that morpholog- 
ical types are not always uniquely assigned (Lahav et al. 
1995), and often one only knows roughly what the actual 
type is, we therefore report the above parameters and ratios 
for the following three disc galaxy classes: lenticular galax- 
ies (SO, so/a); early-type spiral galaxies (Sa, Sab, Sb); and 
late-type spiral galaxies (Scd, Sd, Sm). To help reduce cross 
contamination, the Sbc and Sc types are not used in our 
galaxy class classification. 

In Table|3]one can find the dust-corrected disc and bulge 
magnitudes obtained using equations |3] and |4j respectively. 
While equation |4] was applied to the bulge magnitudes in 
Table 131 equation [3] was not applied to the disc magnitudes 



^ Differences between the Ks filter and the X-band have been 
ignored. 

* Hubble (1926) referred to Sa galaxies as early-type, Sb as in- 
termediate, and Sc as late-type. 
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Figure 3. Standard (if-band) disc and bulge structural proper- 
ties as a function of galaxy Type. For the disc, /lo, h and M^isc a-re 
the central surface brightness, scale-length and absolute magni- 
tude. For the bulge, fic, Ren and Mijuigc are the effective surface 
brightness, effective half-light radius, Sersic index and absolute 
magnitude. The data have been taken from Table [S] except for 
the bulge magnitudes which have been taken from Table |4] For 
each T-Type, the median is marked with a circle and the 'error 
bars' denote the 16 and 84 per cent quartiles of the full distribu- 
tion, rather than uncertainties on the median value. 



from this table but the disc magnitudes obtained using the 
observed disc scalelengths and central surface brightnesses. 
Given that we now have two estimates of the corrected disc 
magnitude, one in Table |3] and the other in Table |4l these 
are compared in Figured and shown to agree within a cou- 
ple of tenths of a magnitude or better. A change of 0.2 mag 
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Table 3. X-band structural parameters of disc galaxies. The median, +/— 34 per cent of the distribution about the median, is shown as 
a function of galaxy type (Column 1). Columns 2 and 3: Disc central surface brightness and scalclongth corrected using equations [T] and 
[21 respectively. Column 4: Disc magnitude computed using equation [3 and the entries in column 2 and 3. Columns 5, 6 and 7: Observed 
bulge effective surface brightness, radius, and Sersic index. Column 8: Bulge magnitude computed using equation |8] and the entries in 
column 5-7 (no dust correction has been applied here). Column 9: Logarithm of the bulge luminosity density pe = pif = Re)^0 k pc~'^. 
Multiplying pe by the appropriate X-band stellar mass-to-light ratio will give the stellar mass density at r = _Rc (see equation [lO} . The 
final three rows are such that T < are the lenticular galaxies, T = 1-3 are the early-type spiral galaxies, and T = 6-9 are the late-type 
spiral galaxies. 
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for the disc magnitude corresponds to a change in log(i3/i3) 
of 0.08. We are also able to show such a comparison for the 
B- and 7-band, as both sets of dust corrections are available 
in these bands. The slight difference in Figure |4] is thought 
to arise from the exponential function not providing a per- 
fect description of dusty discs (see MoUenhofT et al. (2006, 
their Figure 2). Indeed, due to the prevalence of centrally lo- 
cated dust, the best-fitting exponential model can be seen in 
MoUenhofT et al. to overestimate the observed fiux at small 
radii. Therefore, while the corrections between the measured 
and intrinsic central surface brightness and disc scalelength 
given in equations [1] and [2] are appropriate, application of 
equation [T] using /io.intrin from equation [T] and /iintrin from 
equation [2j can lead to a slight over-estimate of the actual 
observed disc magnitude. As a result, the corrections to the 
observed disc magnitude using this approach are not quite 
as large as they should be (C.Popescu & R. Tuffs 2008, priv. 
comm.). 



4-1.1 The bulge-to-disc size ratio Rc/h 

Using exponential bulge models, de Jong (1996b) suggested 
that the bulge-to-disc size ratio was independent of Hubble 
type. This unexpected result was reiterated by Courteau, 
de Jong & Broeils (1996) using an additional ~250 Sb/Sc 
galaxies imaged in the r-band. Adding to the mystery, Gra- 
ham & Prieto (1999) pointed out that de Jong's early- type 
spiral galaxies actually had an average Rs/h ratio — derived 
using his exponential bulge models — that was smaller (at 
the 3(7 level in the K-hand) than the average Re/h ratio 
from his late- type spiral galaxies. 



■24 - 



" -22 



-20 



18 



S 

A 
★ 










m 













• 


A 


• 



• 






A 




▲ 


A 
▲ 


A 
▲ 


8 





▲ 












• 




t 






A 






★ 


★ 


★ 


▲ 




i 










★ 





* 



A 
▲ 



8 







5 

Type 



10 



Figure 4. Comparison of our two estimates of the (corrected) 
disc magnitudes. The filled symbols show the median values ob- 
tained using the corrections in equations [Tl and [2l while the open 
symbols show the median values obtained using equation \3\ B- 
band: stars, /-band: triangles. A'-band: circles. 



Rather than only using the effective radii from de Jong's 
exponential bulge models, Graham & Prieto explored use of 
the Rc values from de Jong's "best-fitting" bulge model; 
which was either an i?i''* model, an R^^^ model or an ex- 
ponential J?i''i model. Doing so, the mean Rc/h ratio was 
shown to be slightly larger for the early-type disc galax- 
ies than the late-type disc galaxies. Following up on this, 
Graham (2001a) fitted Sersic bulge (plus exponential disc) 
models to all of the hght profiles used by de Jong. Graham 



Inclination- and dust- corrected galaxy parameters 9 



Table 4. K band parameters of disc galaxies. The median, ±68/2 per cent of the distribution on either side of the median, is shown as a 
function of galaxy type (Column 1). Column 2: Bulge-to-disc size ratio using columns 6 and 3 from Table [3] Column 3; Bulge magnitude, 
obtained from the observed flux, corrected using equation |4] Column 4: Disc magnitude, obtained from the observed flux, corrected using 
equation [3] Column 5: Bulgc-to-disc flux ratio using columns 2 and 3. Column 6: Number of data points. The bulge-to-total flux ratio 
can be obtained from the expression B/T = [1 + {D/ B)]~^ , where B/D is the bulge-to-disc flux ratio. 
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revealed little difference between the early- and late-type 
disc galaxy bulge-to-disc size ratios, and therefore ultimately 
reached the same conclusion as Courteau et al. (1996) but 
on somewhat different grounds. To explain this apparent dis- 
crepancy with what one sees when looking at spiral galax- 
ies of different type, i.e. that early-type spiral galaxies ap- 
pear to have larger bulge-to-disc size ratios than late-type 
spiral galaxies, Graham (2001a, his figure 21) subsequently 
presented the ice-berg' model in which the intensity of the 
bulge varies (is effectively raised or lowered) relative to the 
intensity of the disc. 

Using ~400 disc galaxies, Table|4]and FigureO present 
the if-band Re/h size ratio as a function of disc galaxy 
type. One of the nice features of this data set is in fact its 
heterogeneous nature. As a result, possible biases in any one 
paper's modelling are minimised and/or effectively cancelled 
by those from the other papers. The results shown here are 
therefore something of a consensus from many papers. One 
can see that the median value is rather stable at around 
0.221q Q3, with the only departure from this being (i) an in- 
crease to 0.31 for the Sa galaxies and (ii) a value of 0.12 for 
T-types less than zero. The first result arises from some of 
the shallow 2MASS data analysed by Dong & De Robertis 
(2006), which is responsible for the scatter to higher Rc/h, 
and also higher log B/D, ratios than displayed by the other 
data. Excluding this data set, the remaining Sa galaxies have 
Re/h = 0.25 with a variance in this distribution of 0.04. 
The second result above is unusual: the size ratio is equal 
to or less than the lower la limit from every other galaxy 
type (not to mention that lenticular galaxies are commonly 
thought to have the largest bulge-to-disc ratios). This re- 
sult is however again solely due to the data from one study. 
Modelling the bulge and bare (and disc) as separate com- 
ponents, Laurikainen et al. (2005) obtained these smaller 
ratios. In contrast, the barless galaxies from Balcells et al. 
(2003) with T-type less than have {R^/h) ~ 0.2. From the 
K-hand panel in Figure [2] one can see that uncertainties in 
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Figure 5. Top panel: X-band, bulge-to-disc size ratio as a func- 
tion of galaxy Type. Bottom panel: Logarithm of the A'-band, 
bulge-to-disc flux ratio as a function of galaxy Type. Values have 
been taken from Tables |4] For each Type, the median value is 
marked with a large circle and the 'error bars' denote the 16 and 
84 per cent quartiles of the distribution; they thus enclose the 
central 68 per cent of each distribution. 



the disc inclination for the Balcells et al. (2003) sample will 
not change this result by much. 

We are unable to include the lenticular galaxies from 
Barway et al. (2007) as neither their fits nor bulge/disc pa- 
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rameters have yet been released. These authors claim that 
the Re/h ratio is not (roughly) constant for all disc galaxies 
and that luminous lenticular galaxies can have values as high 
as 5 to 10, implying the presence of small embedded discs 
whose light peters out well before the half-light radius of the 
bulge component. Such potentially new objects with bulge- 
to-disc size ratios ~50 times greater than the median value 
of regular lenticular galaxies studied by others would surely 
be an important clue to the transition between disc and el- 
liptical galaxies. If real, these discs might be more akin to 
the nuclear discs seen in many elliptical galaxies (e.g. Rest et 
al. 2001). However, we note that the sky-background should 
always be derived independently of one's fitted models. Bar- 
way et al.'s (undesirable) treatment of the sky-background 
as a free parameter when fitting the bulge and disc models 
may be responsible for the extreme size ratios they report. 
Such an approach to modelling galaxies modifies the real 
surface brightness profile, and may result in the occurrence 
of item (vi) in Section [3.11 



4-1.2 The bulge-to-disc flux ratio 

The inclination- and dust-corrected if-band bulge-to-disc 
fiux ratios, B/D, are provided in Table|4]and shown in Fig- 
ure [5] as a function of galaxy type. Unlike the size ratio, the 
fiux ratio does indeed decrease with increasing galaxy type. 
The early-type spiral galaxies (Sa-Sb) have a median value of 
log(_B/_D) = -0.49 (i.e. B/D = 0.32, or equivalently a bulge- 
to-total fiux ratio B /T — 0.24), where B and D represent 
the luminosity (not magnitude) of the bulge and disc respec- 
tively. The late-type spiral galaxies (Scd-Sm) have a median 
value of log(B/D) = -1.40 (B/D = 0.04, B/T = 0.04). 

The tabulated if-band B/D fiux ratios for the 12 lentic- 
ular galaxies in Andredakis et al. (1995. their table 4) give 
a mean B/T ratio of 0.28. Similarly, as reported in Bal- 
cells et al. (2007), the mean (plus or minus the standard 
deviation) of the B/T fiux ratio for the lenticular galax- 
ies analysed in Balcells et al. (2003) is 0.25 (±0.09). Lau- 
rikainen et al. (2005) also obtained a similar result from 
their Tf-band data, reporting a mean value of 0.24 (±0.11) 
from a sample of 14 SO galaxies while the data in Gadotti 
(2008) yields an i?-band value of 0.28 from their 7 lenticular 
galaxies. Our result implies that (50±68/2=) 84 percent of 
lenticular galaxies have B/T ratios smaller than ~l/3. The 
remark in Kormendy (2008) that "almost no pseudobulges 
have B/T > 1/3" could thus be expanded to read "most 
disc galaxies have B/T < 1/3". 

The above ratios, obtained from the collective aver- 
age of nine modern studies (see Table [2]) , are smaller than 
those reported in the 1980s. Indeed, many late-type spiral 
galaxies have bulge fiuxes which are less than 4 per cent of 
their galaxy's total light. While Shields et al. (2008) claim 
that the T=6 (Scd) galaxy NGC 1042 may be a bulgeless 
galaxy, its optical light profile in Boker et al. (2003, their 
Figure 1 jf| together with its JC-band light profile in Knapen 



^ The bump in Bokcr ct al.'s (2003) optical light profile for 
NGC 1042 at 4 ± 2 kpc is due to prominant inner spiral arms 
which cause a clear excess above the disc's underlying exponen- 
tial light profile. 



et al. (2003), plus the diffuse glow of central light in the op- 
tical images, reveals an excess of fiux over the inner 2 kpc 
which is well above that defined by the inward extrapolation 
of the outer exponential light profile. Our corrected if-band 
bulge magnitude (and B/T ratio) for MGC 1042 is -19.13 
mag (0.054), implying a supermassive black hole mass of 
3 X lO^AfQ (Graham 2007), in agreement with the upper 
bound estimated by Shields et al. (2008). We therefore cau- 
tion that it can be difficult when determining if a galaxy is 
actually bulgeless, which is a topic of particular interest at 
the low-mass end of relations involving supermassive black 
hole masses (e.g. Satyapal et al. 2007, 2008). It is not yet 
determined if massive black can holes form (at any redshift) 
before or without a host bulge. Local galaxies may shed in- 
sight into this 'chicken-egg' problem of whether the bulge or 
black hole formed first, if not in tandem. 

Laurikainen et al. (2007, their figure 1) recently pre- 
sented a set of K- and i/-band B/T flux ratios versus galaxy 
type for a sample of lenticular galaxies and the spiral galax- 
ies from the Ohio State University Bright Galaxy Survey 
(Eskridge et al. 2002) . Although that data has not been fully 
corrected for dust attenuation — which is not insignificant 
in the //-band (e.g. Peletier & Willner 1992; Driver et al. 
2008) — the trend which they show embodies the general 
behaviour observed here. In a separate study, Dong & De 
Robertis (2006) claim that their /f-band flux ratio versus 
galaxy type diagram (their fig. 7) agrees with the original 
/3-band results given by Simien & de Vaucouleurs (1986). 
While roughly true, such an agreement should not occur 
given the differing bulge and disc stellar populations (and 
hence colours in the B- and /f-band) for late-type disc galax- 
ies. Aside from this problem, much of the agreement arises 
from the partial cancellation of two significant errors: Simien 
& de Vaucouleurs over-estimated their bulge fiux due to the 
use of i?^'''* models while at the same time they neglected 
the dimming effects of dust. Working with /?^/"-bulge mod- 
els in the /t-band, Dong & De Robertis (2006) overcame 
these issues. 

Recent claims that the luminosity of the bulge, rather 
than the B/D ratio, is the driving force behind the trend 
seen in Figure [5] (e.g. Trujillo et al. 2002; Scannapieco & 
Tissera 2003; Balcells et al. 2007) refiect the results in 
Yoshizawa & Wakamatsu (1975, their figures 1 and 2) and 
echo the remarks in Ostriker (1977) and Meisels & Ostriker 
(1984) that the bulge luminosity, and ergo mass, may be a 
key parameter which distinguishes galax;ies. In Figure |3] we 
show, separately, the bulge and disc magnitude as a func- 
tion of galaxy T-type. One can see that the disc magnitude 
is roughly constant from T= —1 to T= 4 and then it falls by 
a couple of magnitude upon reaching T= 9 (Sdm galaxy). 
On the other hand, the bulge magnitude is roughly con- 
stant from T= —1 to T= 3, but then falls five magnitude 
by T= 9. We therefore confirm that the bulge magnitude 
is indeed predominantly responsible for the trend between 
the bulge-to-disc luminosity and disc galaxy type. In passing 
we note that Hernandez & Cervantes-Sodi (2006) have advo- 
cated that the spin parameter is the physical quantity which 
determines a disc galaxy's Hubble Type (see also Foyle et 
al. 2008); implying that the spin parameter must therefore 
be connected with the magnitude of the bulge. 

Simien & de Vaucouleurs (1986) presented a Figure 
showing the mean difference between the bulge and galaxy 
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magnitude as a function of morphological T-type. Also 
shown in their Figure was the associated statistical uncer- 
tainty on their mean differences, which is of course much 
smaller than the actual scatter about each mean. They fit- 
ted a quadratic relation to this SM-T data, which has been 
popular in semi-analytical studies in which bulge-to-disc ra- 
tios are known but little information about the spiral arms 
is available. Although one could fit a relation to the median 
B/D ratios as a function of T-type in Figure O we feel that 
this may be misleading. For example, if log(i3/_D) = —0.5, 
the actual T-type could range from -1 (SO) to 6 (Scd), while 
an equation would give only one value. The range of values 
reminds us that the bulge-to-disc fiux ratio was never used as 
the primary criteria for classifying disc galaxy morphology. 
However, given the prominace of a bulge reflects the relative 
dominance of certain construction processes, which may be 
different to those which generate spiral patterns, then it may 
at times be desirable to employ such quantitative measures 
without recourse to Hubble type. 

Knowing the stellar mass-to-light ratio of the bulge and 
disc components allows one to transform the bulge-to-total 
flux ratios, B/T, into stellar mass ratios. This aids compar- 
ison with cosmological simulations of galaxies. The stellar 
bulge-to-total mass ratio is given by {l + x[{B /T)~^ — 1]}"^, 
where x = (A//I/)disc/(M/L)buigc is the ratio of the disc and 
bulge stellar mass-to- light ratios, and B/T is of course the 
bulge-to-total flux ratio. If we use a slightly extreme (low) 
value of X ~ 1/2 (see Bruzual & Chariot 2003), then for 
B/T = 0.25 (0.33) we have a bulge-to-total stellar mass 
ratio of 0.4 (0.5). 

To date, most simulations have a tendency to produce 
bulges rather than (pure) discs because of baryon angular 
momentum losses during merger events (Navarro & Benz 
1991; Navarro & White 1994; van den Bosch 2001; D'Onghia 
et al. 2006). To avoid this issue, while still maintaining the 
hierarchical merger trees for CDM halos, potentially unre- 
alistic amounts of feedback have been invoked in the past. 
Early feedback could blow much of the baryons out of high-z 
galaxy discs for their late-time, and necessarily post major- 
merger, re-accretion into some preferred plane and thereby 
produce the disc-dominated galaxies observed in the local 
Universe (see Abadi et al. 2003a). From the ~1000 spiral 
galaxies simulated by Croft et al. (2008), ai z = 1 only 3.5 
per cent have a bulge-to-total stellar mass ratio smaller than 
0.5. Given that 68 per cent of real lenticular galaxies have 
B/T smaller than ~l/3, these simulations are clearly in- 
consistent with what is observed in the Universe. The disc 
galaxy simulated by Abadi (2003b) has a stellar bulge-to- 
total mass ratio of 0.71, and as such it also does not rep- 
resent a typical disc galaxy. While the high mass and force 
resolution ACDM simulated spiral galaxy in Governato et al. 
(2004) has B/T = 0.26, subsequent high- resolution simula- 
tions (e.g. D'Onghia et al. 2006) reveal that resolution issues 
are not the key problem but rather some primary physics 
may still be missing. Current implementation of supernovae 
feedback, for example, has been shown to play a significant 
role in the evolution of simulated galaxies (Okamoto et al. 
2005; Governato et al. 2007; Scannapieco et al. 2008 and 
references therein). 

The current data may hopefully help determine how 
bulges should be assigned their mass in semi-analytical mod- 
els. Such studies build galaxies by a trial and error process 



in which the model assumptions, parameters and processes 
are collectively tweaked until the simulations resemble real 
galaxies. For example, when a simulated disc becomes un- 
stable (Mo et al. 1998; Cole et al. 2000), Croton et al. (2006, 
their equation 22) transfer enough disc mass to the bulge in 
order to maintain stability, while Bower et al. (2006, their 
equation 1) transform the entire system into an elliptical 
galaxy. In the case of minor mergers (mass ratio <0.3), some 
works assign all the accreted gas to the disc while transfer- 
ring all the accreted stars to the bulge (de Lucia et al. 2006), 
while others add both the stars and cold gas to the disc of 
the primary galaxy (Kauffmann & Haehnelt 2000). The new 
bulge-to-disc flux ratio data may help to reflne how these 
semi-analytical models operate. 



4-1.3 Bulge profile shape 

Although equations to correct the Sersic parameters pertain- 
ing to the bulge are not yet available, we note that the influ- 
ence of dust at wavelengths around 2 microns is noticeably 
less than in the optical bands (e.g. Whitford 1958; Cardelli, 
Clayton & Mathis 1989; Fitzpatrick 1999). The signiflcant 
change to ^,o,k with inclination which is seen in Figure [5] is 
not because of dust but predominantly due to the greater 
line-of-sight depth through the (near-transparent) disc with 
increasing disc inclination. In the absence of dust, a spheri- 
cal bulge will look the same from all angles, i.e. inclinations. 
We therefore do not expect the if-band bulge parameters to 
change greatly (i.e. more than a few tenths of a magnitude, 
or ~20 per cent in size) with inclination, and we report their 
observed values in Table [3] 

As can be seen in Figure[31 Sc galaxies and earlier types 
typically have values of n ^ 2 and sometimes as high as 4. 
While Scd galaxies and later types tend to have values of 
n < 2, some Sd galaxies can still have values of n greater 
than 2. It may therefore be misleading to label all Sc and 
later-type galaxies as possessing pseudobulges built from ex- 
ponential discs (Erwin et al. 2003; Kormendy & Kennicutt 
2004, their section 4.2; Athanassoula 2008). While a value 
of n = 1 in Sersic's model describes an exponential profile, 
it does necessitate the presence of a fiattened, rotating, ex- 
ponential disc rather than a spheroidal, pressure supported 
bulge. Dwarf elliptical galaxies also have values of n < 2 (e.g. 
Caon et al. 1993; Young & Currie 1994; Jerjen et al. 2000) 
— they define the low-mass end of a linear relation between 
magnitude and Sersic index (e.g. Graham et al. 2006b, their 
figure 1; Nipoti, Londrillo & Ciotti 2006). Dwarf ellipticals, 
however, are generally not considered to have formed from 
the secular evolution of a disc0 Subsequently, disc galaxy 
bulges with values of n < 2 (or n ~ 1) are not necessar- 
ily pseudobulges and should not be (re-)classified as such 
based on this criteria alone. To clarify Kormendy & Ken- 
nicutt (2004), while pseudobulges may have n ~ 1, a bulge 
with 71 ~ 1 need not be a pseudobulge. 

Although not a new result, the preponderance of bulges 
with n < 4 combined with the prevalence of late-type disc 



^" Note: even if dwarf elliptical galaxies had formed from the 
violent redistribution of a disc, they need not preserve the disc's 
exponential radial profile shape. 
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galaxies having bulges with n ~ 1 necessitates that simu- 
lations of disc galaxies neither use nor create an i?^'^*-like 
bulge for every galaxy. While the projection of Hernquist's 
(1990) useful model — employed in many numerical simu- 
lations (e.g. Springel et al. 2005b; Di Matteo et al. 2005) 
— reproduces an R^^'^ profile, it does not provide a particu- 
larly good representation of bulges, such as the Milky Way's 
bulge, which have n ~ 1 (e.g. Terzic & Graham 2005, their 
Figure 7; Terzic & Sprague 2007). Furthermore, efforts to 
simulate a range of disc galaxies with different morphologi- 
cal types by only varying the bulge-to-disc mass ratio, but 
using a Hernquist model for every bulge (e.g. Springel et 
al. 2005a), i.e. effectively assigning a de Vaucouleurs profile 
to every bulge, will not generate a realistic set of galaxies. 
In the case of bulges formed from disc instabilities, one re- 
quires the production of bulges that have exponential light 
profiles rather than Hernquist-like density profiles. In regard 
to simulations of minor mergers, in which bulges are built up 
within discs, the observed correlation between bulge profile 
shape and luminosity, and also supermassive black hole mass 
(Graham & Driver 2007), must be reproduced. The recov- 
ery of bulge parameters using i?^'''*-based models is known 
to introduce systematic biases with profile shape and thus 
also black hole mass (e.g. Trujillo et al. 2001, their Fig- 
ures 4 and 5; Brown et al. 2003, their Figure 7). Therefore, 
once realistic bulges are generated, the computation of their 
sizes and binding energy by assuming a Hernquist model 
would be inappropriate and lead to spurious correlations 
when constructing, for example, a black hole fundamental 
plane (Younger et al. 2008) . 

Density models such as that from Einasto (1965) and 
Prugniel & Simien (1997) contain a "shape parameter" that 
effectively captures the range of structural shapes (Sersic 
indices) which real bulges are observed to possess Such 
models have already provided insight into the gravitational 
torques which oblate and triaxial bulges generate. Interest- 
ingly, it is bulges with smaller Sersic indices which have 
a greater non-axisymmetric gravitational field: this is be- 
cause bulges increasingly resemble central point sources as 
the Sersic index increases (Trujillo et al. 2002). 

By construction, the three-parameter (pe, Rc,n) density 
model of Prugniel & Simien (1997; their equation B6) con- 
tains two identical parameters to Sersic's model and can be 
written as 



p(0 = P.(^)"^-^[<^^^"'^^"^, 

_ M h 6"(i-f' r(2n) 
^"^T Wc 2r(n(3-p)) ' 



(9) 
(10) 



where r is the internal, i.e. not projected, radius (see 
Marquez et al. 2001; Terzic & Graham 2005, their equa- 
tion 4). For clarity, the subscript n has been dropped from 
the term 6„. The new parameter is the internal den- 
sity at r = Tie and provides the normalisation such that 
the total mass from equation ((9} equals that obtained from 



Kincmatical expressions associated with these models, plus 
equations to convert projected densities such as {fi)c into three- 
dimensional densities, are given in Terzic & Graham (2005), Mer- 
ritt et al. (2006) and Graham et al. (2006a). 



Table 6. V band parameters. The median, ±68/2 per cent of the 
distribution on either side of the median, is shown as a function of 
galaxy type (Column 1). Columns 2 and 3: Face-on, disc central 
surface brightness and scalelength corrected using equations [T] 
and [21 respectively. Column 4: Disc magnitude computed using 
equation [7] and the entries in column 2 and 3. Column 5: Number 
of data points. 
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equation ([5} after applying the appropriate stellar mass-to- 
luminosity ratio M/L. Equation (|10p has the same func- 
tional form as Sersic's model although multiplied by an adi- 
tional power- law term with exponent p = p{n). We adopt 
Lima Neto et al.'s (1999) expression for p, for which a high- 
quality match between deprojected Sersic profiles and the 
above expression is obtained when p = 1.0 — 0.6097/n -I- 
0.05563/n2. 

After converting the observed A'-band bulge effective 
radii into units of parsecs, and the observed K-hand effective 
surface brightnesses into units of absolute solar luminosity 
per square parsec, we have used equation (|10|l to derive the 
mean pc/{M/L) values for bulges of each galaxy type and 
class. These luminosity densities are provided in Table [31 
To do this we used an absolute iC'-band magnitude for the 
Sun of 3.33 mag (Cox 2000). Multiplying by one's preferred 
stellar M/L ratio gives the stellar mass density of the bulge 
at r = Re. The mean Rc, n and pc values can be used to 
construct realistic bulges using the Prugniel-Simien model, 
and/or to test if simulated bulges match the typical density 
profiles of real galaxy bulges. 



4.2 The optical 

As with the TiT-band data, we have collated literature data 
on the structural parameters of galaxies measured at opti- 
cal wavelengths (see Table [2}. After first removing any in- 
clination corrections applied in these papers, we again em- 
ployed equations [T] and [2] to determine the disc central sur- 
face brightnesses and scalelengths, and we computed the 
disc magnitudes using equation [T] The results are shown in 
Tables [5] to H 

While dust is expected to modify the individual bulge 
Sersic parameters at optical wavelengths — and for this rea- 
son we don't present such values here — we can apply the 
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Table 5. B band parameters. The median, ±68/2 per cent of the distribution on either side of the median, is shown as a function of 
galaxy type (Column 1). Columns 2 and 3: Face-on, disc central surface brightness and scalelength corrected using equations [Tl and [2l 
respectively. Column 4: Disc magnitude computed using equation [7] and the entries in columns 2 and 3. Column 5: Bulge magnitude, 
obtained from the observed flux, corrected using equation |4] Column 6: Disc magnitude, obtained from the observed flux, corrected using 
equation O Column 7: Bulge-to-disc flux ratio using columns 5 and 6. Column 8: Number of data points. 
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Table 8. I band parameters. The median, ±68/2 per cent of the distribution on cither side of the median, is shown as a function of 
galaxy type (Column 1). Columns 2 and 3: Face-on, disc central surface brightness and scalelength corrected using equations [Tl and l2l 
respectively. Column 4: Disc magnitude computed using equation [7] and the entries in columns 2 and 3. Column 5: Bulge magnitude, 
obtained from the observed flux, corrected using equation |4] Column 6: Disc magnitude, obtained from the observed flux, corrected using 
equation [3] Column 7: Bulge-to-disc flux ratio using columns 5 and 6. Column 8: Number of data points. 
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magnitude corrections from equation|4]to the observed bulge 
magnitudes. Table [5] sho-ws these B-band bulge (and disc) 
magnitudes and the logarithm of the ratio of their luminosi- 
ties. This latter quantity is plotted in Figure [6] The B-band 
flux ratios obtained from Simien & de Vaucouleurs (1986), 
which are still used by many research groups (e.g. Springel 
et al. 2001; Mathis et al. 2002; De Lucia et al. 2006) to build 
galaxies and assign galaxy types, are up to 0.5 dex larger. If 
the i?i'''i-derived flux ratios from Simien & Vaucouleurs -were 
corrected for dust, then they would be greater still. While 
our galaxy parameters have been corrected for dust, they 



have also been derived using i?i''" rather than R^''^ bulge 
models which, collectively, results in smaller flux ratios. 

Having corrected the B-band B/D ratios for dust does 
not imply that they should be equal to the iC'-band values. 
In fact, differences are expected and will reflect the varying 
stellar populations of the bulges and discs. For example, 
a redder intrinsic (i.e. before the effects of dust) SED for 
a bulge will result in a larger B/D ratio in the 7^-band 
than in the B-band. In Figure [7] we have co-plotted the K- 
and B-band B/D values as a function of galaxy type. The 
SO-Sab galaxies have roughly the same B/D ratio, possibly 
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Table 7. R band parameters. The median, ±68/2 per cent of the 
distribution on cither side of the median, is shown as a function of 
galaxy type (Column 1). Column 2: Bulge magnitude, obtained 
from the observed flux corrected using equation |4] Column 3: 
Disc magnitude, obtained from the observed flux corrected using 
equation \3\ Column 4: Bulge-to-disc flux ratio using columns 2 
and 3. Column 5: Number of data points. 
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Figure 6. Logarithm of the B-band, bulge-to-disc flux ratio as a 
function of galaxy Type (see Table [5] column 7). The dashed line 
traces the B-band values from Simien & de Vaucoulcurs (1986). 



indicating similar stellar population, i.e. if the bulge is red, 
then the disc will also be red (Peletier & Balcells 1996). 
Although we do caution that other factors may be at play. 
Laurikainen et al. (2005) were careful to separate the bar flux 
from the bulge flux in their if-band analysis of early-type 
disc galaxies while Balcells et al.'s (2003) if-band sample 
did not possess bars. The optical studies might have overly 
large BjD ratios due to the assignment of bar flux to the 
bulge. 

The later type spiral galaxies tend to have a B-band 
\og{B/D) value which is nearly some ~0.3 dex smaller than 
the A'-band value. That is, for the intermediate and late- 
type spiral galaxies, the B-band B/D ratio is a factor of 2 
smaller than the A'-band ratio. The obvious interpretation 
is that the later types have bluer — younger and/or more 
metal poor — discs relative to their bulges. Regarding the 
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Figure 7. Median logarithm of the bulgc-to-disc flux ratio as a 
function of galaxy Type. The number of galaxies involved in each 
data point can bo seen in Tables |4] to |8] 



Sm and irregular galaxies, the situation becomes somewhat 
mixed, and the diagram also suffers from low number statis- 
tics at this end. To provide a rough idea of the uncertainty 
on the plotted data points in Figure [T] we note that for a 
distribution with sample standard deviation equal to a, the 
uncertainty on the mean is a / yJ\N), where is the num- 
ber of data points. From Table|3]we have a ~ 0.4, while the 
value of N can be seen in Tables |4] to [8] 



5 BIVARIATE DISTRIBUTIONS 

Scaling relations, if found to exist, can (i) provide insight 
into the physical mechanisms which shape galaxies and 
(a) supply helpful restrictions for simulations which try 
to match the real Universe and thereby understand it. In 
this section we present dust-corrected brightness-size dis- 
tributions for galaxy discs (and also bulges and ellipti- 
cal galaxies). Because the majority of disc galaxies consist 
of two main, physically distinct, components — a three- 
dimensional bulge and a relatively flat, two-dimensional disc 
— we do not treat the disc galaxies as single entities. Instead, 
we separately show bivariate distributions for the structural 
properties of the bulges and the discs. 

The data from Dong & De Robertis (2006) is not in- 
cluded from here on. Their challengingly-shallow (eight sec- 
ond exposure) 2MASS images may be responsible for a num- 
ber of apparent outliers which would slightly bias the fits in 
the following section if they had been included. 



5.1 The ^o-log/i distribution 

Figure shows the Jf-band distribution of face-on, central 
disc surface brightness jio versus disc scalelength h. Both 
of these parameters have been fully corrected for the influ- 
ence of dust. For exponential discs, one can easily transform 
this diagram so that it involves the effective radii, Rc, and 
effective surface brightnesses, pe- This is achieved by using 
^^o.disc = 1.678/1 and /ie.disc = A*o + 1-822, or the mean ef- 
fective surface brightness within i?e,disc, which is such that 

(M)o,disc = /iO + 1.123. 
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0.1 1 10 

Figure 8. X-band. Top panel: Central disc surface brightness 
(corrected using equation [TJ versus disc scalelength h (corrected 
using equation [2)1 . Bottom panel: Bulge effective surface brigfit- 
ness (observed) versus the bulge effective radius Re (observed). 
The upper bright envelope is traced here with the (empirical) 
lines /iQ^bright = 14.85 + 2.2 log h (panel a) and /.tc, bright = 
15.3 + S.Ologiie (panel b). Galaxy types Sbc (T = 4) and earlier 
are denoted by the circles, while later galaxy types are denoted 
by the crosses. 



Due to selection effects, small faint galaxies are likely to 
be missed and thus the lower-left of Figure [8^ may be miss- 
ing galaxies. Tfie upper envelope in Figure [8^ does however 
reflect a real cutofT in the distribution of galactic parame- 
ters, simply because bigger and brighter discs would be seen 
if they existed. The late- type disc galaxies can be seen to oc- 
cupy a much broader region of this diagram than the early- 
type disc galaxies, and as such they are not particularly 
well represented by a linear surface brightness-size relation. 
Moreover, given the diflfering distribution of the early- and 
late- type disc galaxies (see also Figure [9]), one can see that 
it is predominantly the late-type disc galaxies for which we 
may have an incomplete distribution in the lower-left quad- 
rant. 

Figure [9] shows the B-band /^o-log h relation in which 
the same general behaviour is observed. One of the immedi- 
ate results that can be drawn from these figures is that the 
Freeman (1970) law, pertaining to a constant central disc 
surface brightness, requires modification. Not only do low 
surface brightness galaxies exist with fainter surface bright- 
nesses — which has long been known (e.g. Disney 1976, 
1998) — but the bright limit does not simply have the canon- 
ical Freeman surface brightness of 21.65 B-mag arcsec"^. 
Instead, the bright limit is observed to depend strongly on 
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Figure 9. B-band. Central disc surface brightness (corrected 
using equation [T} versus disc scalelength (corrected using equa- 
tion [2]l. Rather than a Freeman (1970) law of constant sur- 
face brightness, an upper bright envelope, which has been seen 
before (e.g. Graham 2001b and references therein), is evident. 
The empirically-fitted dotted line in this diagram is such that 
IJ-o,brighf = 18.6 -1-2.5 log /i. Galaxy types Sbc (T = 4) and earlier 
are denoted by circles, while later galaxy types are denoted by 
crosses. The noticeably broader distribution for the late-type disc 
galaxies has also been noted before (Graham & de Blok 2001). 



disc size such that Sebright ~ 2.5 log h (Graham 2001b; Gra- 
ham & de Blok 2001, and references therein). Obviously, this 
relation (quantified in section [5?3|l also implies an improve- 
ment over the simple use of average disc scalelengths and 
central surface brightnesses to represent each of the (early- 
type) disc galaxy types. 



5.2 The /ic-log i?c distribution 

The 7^-band distribution of bulge effective surface bright- 
nesses /ic and effective radii Rc are plotted in Figure [S]d. As 
with the disc parameters, the bright envelope to this distri- 
bution is real and not an artifact of galaxy selection criteria. 
The bright boundary to the distribution is roughly described 
by the empirical relation ^to, bright ~ 3 log Rc - Figure |8Jd qual- 
itatively agrees with the (surface brightness)- size diagram in 
Binney & Merrifield (1998, their Fig. 4. 52) which was based 
on optical images of 66 galaxies from Kent (1985). The ad- 
vances here are (i) sample size (ii) the influence of dust is 
less in the if-band and (iii) 7?^^" rather than R^^* models 
have been applied to every bulge. 



5.3 The luminosity-size distribution 

Given that L — 2TvIoh^ for exponential discs, the size- 
luminosity diagram is simply an alternative representa- 
tion of the more traditionally used (surface brightness)- 
size diagrams seen in section [STTI For convenience, both are 
shown here, although we emphasize that the upcoming size- 
luminosity relations could have been equally as well obtained 
by fitting lines to the data in Figures and[^ The dashed 
lines in those figures have been mapped into Figure [10] which 
shows the (dust-corrected) disc magnitudes versus the (dust- 
corrected) disc scalelengths. 

The early-type disc galaxy distribution in Figure lTOl has 
been fit with a linear relation such that 
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log/i = a- ^(Mdisc + Afo), (11) 

and thus h oc L''. The quantity Mo is simply a constant to 
help centralise the fit and thereby reduce the uncertainty 
on the intercept a. We have employed Feigelson & Babu's 
(1992) code SLOPES to perform six regression analyses: 
ordinary least squares of Y on X, OLS(Y|X); and the in- 
verse OLS(X|Y); the line which bisects these two; minimisa- 
tion of the perpendicular residuals (known as orthogonal re- 
gression); the geometric mean of OLS(Y|X) and OLS(X|Y), 
referred to as the reduced major axis regression line; and 
finally the arithmetic mean of the two OLS lines. Here, 
Y — log/i and X = Afdisc, and the results of the regres- 
sion are shown in Tables l9l and 1101 Table [TOl also shows the 
results of the regressions for different galaxy types, although 
only in the /("-band due to the numbers of galaxies involved. 

From Graham & Driver (2005, their equation 12) one 
has that the absolute magnitude 

M = - 2.5 log[27ri?^] - 36.57, (12) 

with Rc in kpc. Therefore, for a disc with n = 1, 

Mdisc = (/io + 1.123) - 2.5 log[27r(1.678/i)^] - 36.57. (13) 

From this relation one can readily express the disc size- 
luminosity relations (Tables |9] and llOp in terms of size and 
surface brightness. For example, combining eguations llll and 
ll3l to eliminate Mdisc, one has 

56 — 2.5, , 2.5a 
MO = ; log h - Mo + —— + 38.57. (14) 



If the typical central surface brightness = —2.5 log Io 
was constant for the discs, although we know from Figures [8] 
and [9] that it is not, then one would find h ~ While 
a symmetrical regression of the optical data suggests h oc 
^o.75±o.04 ^jjg ^jg^g ^£ early-type disc galaxies (Tabled, 
the near-infrared data has a slope which is smaller, albeit 
only at the ~1.5a level. 

A quick visual inspection of Figure [10] reveals that fit- 
ting to the full data set, i.e. both the early- and late-type 
disc galaxies, would produce h ^ relations with shallower 
slopes (smaller b values) than those reported in Table [9] We 
have not done this because (i) sample selection effects may 
have introduced an artificial boundary to the distribution 
of the late-type disc galaxies and (ii) from the — log h 
diagrams it was already clear that no real relation exists for 
the late-type disc galaxy population. 

Before contrasting our results with others, we note 
that qualitative differences should not be surprising. We 
have derived intrinsic, fully dust corrected, disc scalelengths 
and luminosities from Sersic i?^''" bulge plus exponential 
disc decompositions. We first compare our K-hand results 
with Courteau et al. (2007) who presented isT-band galaxy 
luminosity-size relations for 360 2MASS galaxies binned into 
different galaxy types. At least for systems with small bulge- 
to-disc ffux ratios, the distribution of disc light is roughly 
approximated by the distribution of the galaxy light. In- 
deed, as shown in Graham (2002, his Figure 6), galaxy half 
light radii are roughly equal to the half-light radii of the disc 
component. 

Courteau et al. (2007) used an orthogonal regression 
analysis which can be seen in Tables [9] and [10] to give a 
shallower relation to our data than the three other symmet- 
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Figure 10. Inclination- and dust-corroctcd magnitude versus 
disc scalclcngth (via equation [2J . Galaxy types Sbc (T = 4) and 
earlier are denoted by the circles, while later galaxy types are 
denoted by the crosses. The differing distributions between the 
galaxy types in these figures is directly related to their different 
distributions in the /xo-log/i plane (Figures |8] and |9} . Similarly, 
evidence of a relation here is simply a re-representation of evi- 
dence for a relation in the fiQ-log h plane. The fitted lines to the 
discs of the early-type disc galaxies are given by the bisector re- 
gression of equation 1111 (see Tabic The dotted lines seen here 
denote the bright boundaries shown in Figures |8ji and |9] 
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Table 9. Linear regression coefficients for the size-luminosity relations of disc galaxies with T < 4 (i.e., <Sbc). The fitted relation is 
such that log/i = a — 0.4b(A'/(jisc + Afo). The OLS bisector lines of regression are shown in Figure [TOl 

Method B-band (Mq = 20) V-band (Mq = 21) /-band (Mq = 22) K-band (Mo = 24) 

abababab 

OLS(Y|X) 0.40 ± 0.02 0.65 ± 0.04 0.48 ± 0.03 0.59 ± 0.08 0.49 ± 0.02 0.65 ± 0.05 0.54 ± 0.01 0.49 ± 0.04 

OLS(X|Y) 0.34 ±0.03 0.86 ± 0.06 0.42 ± 0.04 1.00 ±0.14 0.46 ± 0.03 0.89 ± 0.07 0.54 ± 0.02 0.93 ± 0.07 

OLS bisector 0.37 ± 0.02 0.75 ± 0.05 0.46 ± 0.03 0.77 ±0.09 0.47 ± 0.02 0.76 ± 0.05 0.54 ± 0.02 0.69 ± 0.03 

Orthogonal 0.38 ± 0.02 0.71 ± 0.05 0.47 ± 0.04 0.71 ±0.11 0.48 ± 0.02 0.73 ± 0.06 0.54 ± 0.01 0.59 ± 0.04 

Reduced major axis 0.37 ± 0.02 0.75 ± 0.05 0.46 ± 0.03 0.77 ±0.09 0.47 ± 0.02 0.76 ± 0.05 0.54 ± 0.02 0.67 ± 0.04 

Mean OLS 0.37 ± 0.02 0.75 ± 0.05 0.45 ± 0.03 0.79 ± 0.09 0.47 ± 0.02 0.77 ±0.05 0.54 ± 0.02 0.71 ± 0.04 



Table 10. Linear regression coefficients for the A'-band disc size-luminosity relations of different disc galaxy types: log h = a — 
0.46(Mdi,c + 24). 
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0.61 


± 


0.01 


0.52 


± 


0.03 


0.73 


± 


0.03 


0.40 


± 


0.04 


Orthogonal 


0.51 


± 


0.03 


0.70 


± 


0.09 


0.55 


± 


0.02 


0.58 


± 


0.05 


0.58 


± 


0.02 


0.39 


± 


0.04 


0.62 


± 


0.03 


0.23 


± 


0.05 


Reduced major axis 


0.51 


± 


0.03 


0.78 


± 


0.08 


0.54 


± 


0.02 


0.67 


± 


0.04 


0.60 


± 


0.01 


0.50 


± 


0.03 


0.70 


± 


0.03 


0.36 


± 


0.03 


Mean OLS 


0.61 


± 


0.03 


0.84 


± 


0.09 


0.54 


± 


0.02 


0.70 


± 


0.04 


0.61 


± 


0.01 


0.53 


± 


0.04 


0.74 


± 


0.04 


0.42 


± 


0.05 



rical regression analyses which we used. For our (potentially 
sample-selection biased) late-type galaxies (T > 5) we have 
the 7^-band relation h ~ j^o.23±o.05 (orthogonal regression), 
which is in fair agreement with Courteau et al.'s relation 
igaf "^" ^ for Sd galaxies (their Table 3). For our Sc, Sb 
and Sa galaxy samples we find exponents of 0.39 ± 0.04, 
0.58 ± 0.05 and 0.70 ± 0.09 respectively. Such an increase 
is in qualitative agreement with Courteau et al. who report 
values of 0.39 ± 0.03, 0.43 ± 0.04 and 0.57 ± 0.04. 

Courteau et al. additionally presented a size-luminosity 
diagram for 1300 disc galaxies with /-band imaging, report- 
ing h ~ L°-33±o.02 galaxies), h ~ lO-37±o.02 galaxies) 
and h ~ /,0-55±0 05 (Sa galaxies). When we fit aU of our /- 
band data, i.e. late- and early-type disc galaxies, we obtain 
h ~ £^o.45±o.05 (orthogonal regression), while our fit to the 
SO-Sbc galaxies is notably steeper, with h ~ 7;,o.73±o.06 
thogonal regression). Due to our smaller number of galax- 
ies with /-band data, we do not attempt a finer division of 
galaxy type. 

We have also compared our result with the disc galaxy 
size-luminosity relations in Shen et al. (2003) . In their work 
they show the r' galaxy luminosity versus the galaxy half 
light radius. At the luminous end of Shen et al.'s size- 
luminosity diagram (M^/ < —20 mag), they report that 
radius varies with L°'^^. They performed an OLS(Y|X) re- 
gression analysis on binned data, making it additionally 
difficult to perform an exact comparison with our results. 
At magnitudes brighter than M^i — —20 mag, their Fig- 
ure 5 reveals that the early-type disc galaxies follow the re- 
lation h ~ L"'^*. At lower-luminosities they report shallower 
slopes. 

Shen et al. (2003) assumed that galaxies with concen- 
tration parameters less than 2.86 are disc galaxies. When 
dealing with bright galaxy samples this is a reasonable ap- 
proximation which has been illustrated with visually classi- 
fied samples. However, a problem occurs when one expands 
one's sample to fainter magnitudes because most dwarf el- 
liptical galaxies have n < 2 and thus concentrations less 
than ~ 2.86 (Graham et al. 2005, their table 1). As a re- 



sult, faintward of ~ —19 r'-mag (~ —20 z'-mag) in Shen 
et al.'s Figures 4-7 (and 10), one does not know if the sys- 
tems are actually early- or late-type galaxies. This can how- 
ever be resolved by looking at Shen et al.'s Figures 8-9, in 
which it becomes apparent that the curved nature of their 
disc galaxy size-luminosity relation is due to the mixing of 
galax;y types at the low-luminosity end. This implies that (i) 
the curved size-luminosity relations for late-type galaxies in 
Table 1 from Shen et al. are not appropriate and (ii), the lin- 
ear size-luminosity relations for early-type galaxies in their 
Table 1 are also not appropriate at faint magnitudes. This 
has prompted us to revisit the size-luminosity diagram for 
elliptical galaxies, to which we are able to add and compare 
the distribution of bulges. 

5.3.1 Elliptical galaxies 

From —18 < AIb < —13 mag, the effective radii of dwarf 
elliptical galaxies display no real trend with magnitude but 
instead scatter around a value of ~1 kpc (Smith Castelli 
et al. 2008; Dabringhausen et al. 2008; Figure [TTK). This 
often over-looked result was previously shown in Binggeli & 
Jerjen (1998) and actually also noted in Shen et al. (2003, 
their Figures 8 and 9) but the implications of which were 
not propagated into all of their diagrams. 

The curved distribution of data points in Figure [TTK can 
be predicted (see the line in this figure) from the following 
linear Al — /lo and M — log n (Figure lllb ) relations for el- 
liptical galaxies. Using equations 7, 9 and 12 from Graham 
& Driver (2005), one has an expression involving M,Ra,fJ,o 
and n for any Sersic profile. For elliptical galaxies, the fj,o 
term in this (not shown) expression can be replaced with M 
via the empirical relation 

Mb = ^/io.s - 29.5 (15) 



from Graham & Guzman (2003), to give the relation 



log/?e[kpc] = 1.137+0.2176„ 



Mb 
' 10 



+0.5 log 



nr(2?i)e'' 



,(16) 
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where 6„ ^ 1.9992n - 0.3271, for 0.5 < n < 10 (Capaccioli 
1989). The value of n can also be expressed in terms of Mb 
using the empirical relation 

n = io-(^'^-^+^'b)/9a (^7) 

from Graham & Guzman (2003). This is how the line in 
Figure lllb , was derived. It can be seen to match the be- 
haviour of the data rather well, which has (only) an ap- 
proximately linear behaviour {Rc oc L"'^ from the data in 
Figure [TTk ) at luminosities brighter than Mb ~ — 19 mag. 
The gradual curvature in the bright end of this L — i? re- 
lation is independently supported by the data in Liu et al. 
(2008, their Fig.l2), for which their brightest cluster galaxies 
(when using 25 r-mag/arcsec^ isophotal magnitudes) have 
R50 oc L°'^^ . The exponential-like envelopes of cD galaxies 
(Seigar et al. 2007) does however complicate the derivation 
of reliable brightest cluster galaxy parameters and we do not 
explore this issue here. 

The curved size-luminosity relation for elliptical galax- 
ies, predicted from two linear relations, is probably not a 
fundamental relation. Due to the continually varying slope of 
the L — R relation, the mean slope for dwarf and bright ellip- 
tical galaxies will be different. This is, however, not evidence 
for two disjoint species of galaxies. Similarly, as explained 
in Graham & Guzman (2003), the continually curved size- 
(effective surface brightness) relation for elliptical galaxies 
is not evidence for two distinct types of galaxy. Instead, the 
curved relation is expected from the linear relations which 
dwarf and luminous elliptical galaxies are observed to co- 
define. 

Having derived the curved size-luminosity relation for 
elliptical galaxies from two linear relations, we examine why 
the bright end slope is different to that reported by Shen et 
al. (2003). Shen et al. used Petrosian quantities for the lumi- 
nous elliptical galaxies, however Petrosian magnitudes miss 
the outer flux in galaxies, and increasingly so for brighter 
galaxies. While the Petrosian magnitude is therefore fainter 
than the total magnitude, it is actually the Petrosian half 
flux radii R50 which is most heavily impacted in compari- 
son with the effective half light radii _Ro. This is quantified 
in Graham et al. (2005). Looking at Shen et al.'s Figure 4, 
one can see that their early-type galaxies display a linear 
size-luminosity relation from —19.75 > M^/ > —23.75. From 
Fukugita et al. (1995) we have B — r' — 1.32 and so the above 
magnitude range corresponds to —18.43 > Mb > —22.43. 
From Figure lllb we see that this magnitude range cor- 
responds to a range in Sersic index from 3 to 7. We can 
now use the relation Rso fs [1 - 6.0 x W'^ {R90 / R5of -^^]Re 
from Graham et al. (2005) to quantify the difference be- 
tween the Petrosian half light radii -R50 (used by Shen et 
al.) and the effective half light radii Re (used in Figure [TTk ). 
When n = 7 {Rgo/Rso = 3.62), R50 = 0.42i?c. When n = 3 
{R90/R50 = 3.17), R50 = 0.82_Ro. Such a systematic reduc- 
tion to the radii, as the elliptical galaxy luminosity increases, 
results in a reduction (of a factor log[4.2/0.82] « 0.71) to the 
exponent b in the relation R (x L^, and likely explains why 
Shen et al. (2003) find -R50 <x I/pct°osian- The Sersic radii 
from Blanton et al. (2005) — which Shen et al. additionally 
used — are also known to systematically underestimate the 
true effective half light radii as the Sersic index increases (see 
Figure 9 from Blanton et al. 2005). Underestimation also oc- 
curs when using R^'''^ models to describe systems with n > 4 
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Figure 11. B-band. Panel a) Luminosity-size diagram for el- 
liptical galaxies. The data have come from the compilation in 
Graham & Guzman (2003, see their Figure 9). The line is not a 
fit to the data, but a prediction based upon the linear relations 
between absolute magnitude and central surface brightness (Gra- 
ham & Guzman 2003, their Figure 9f) and absolute magnitude 
and the logarithm of the Sersic index — shown here in panel b) 
for convenience but taken from Graham & Guzman (2003, their 
Figure 10). A fixed change in Sersic index n, for example &n = 1, 
causes a large change to the profile when n is small, but only 
a small change when n is large. Because of this, plots of galaxy 
parameters such as magnitude or colour will appear curved (flat- 
tening at large n) rather than linear, if they are plotted against 
n rather than log n. 

(see Trujillo et al. 2001, their Figure 4) which may explain 
why Bernardi et al. (2003) report that R ~ LP'^^ . There- 
fore, while the above results are not necessarily wrong, it 
does highlight that if we are to compare theory (e.g. Mayer, 
Governato & Kaufmann 2008) and observations, it is impor- 
tant that we measure radii and luminosities in a physically 
appropriate and consistent manner. 

5.3.2 Bulges 

Finally, with new data in hand, the K-hanA size-luminosity 
distribution for the bulges is shown in Figure [T^] No real 
trend is evident and we refrain from fitting a relation to 
this scatter diagram. We can however explore if the upper 
boundary in this Figure is real. Potentially, just as low sur- 
face brightness disc galaxies avoided detection due to ef- 
fectively being hidden by the sky background flux (Disney 
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Figure 12. ii'-band. Size-luminosity diagram for bulges. Galaxy 
types She (T = 4) and earlier are denoted by the circles, while 
later galaxy types arc denoted by the crosses. 
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Figure 13. Size-luminosity diagram for elliptical galaxies (open 
circles) and bulges (dots). The ii'-band bulge magnitudes have 
been shifted to the B-band using a constant B — K colour of 4.0. 



1976; Longmore et al. 1982; Davies et al. 1988), a population 
of large, low surface brightness bulges (which could populate 
the upper left portion of Figure I12|) may be submerged be- 
neath the flux of their associated discs. The existence of such 
bulges is likely to be rare for the simple fact that galaxies in 
which no bulge is detected, i.e. bulgeless galajcies, are rare. 
If a bulge was to exist with, for example, fJe = 3 kpc and an 
absolute ii'-band magnitude of —20 mag, then it would have 
a mean effective surface brightness of 20.95 K-ra&g arcsec^^ 
(Graham & Driver 2005, their eq.l2). If this bulge had a 
Sersic index n = 1, then it would have a Jf-band central 
surface brightness of 19.83 mag arcsec"^. If bright enough 
relative to the host disc, this will result in a central protru- 
sion of the galaxy light profile, above that of the disc profile. 
For late-type, bulgeless disc galaxies (say Sd or later), we can 
see from Table [3] that their central disc surface brightnesses 
are not bright enough to hide such a bulge. On the other 
hand, a bulge with _Re = 3 kpc and an absolute JsT-band 
magnitude of —18 mag, and thus a central surface bright- 
ness of 21.83 mag arcsec"'^, could be hidden. Although we 
note that if the ratio <Rc/h>~ 0.2 — 0.25 holds, it would 
require the late-type disc galaxy has a scalelength around 
12-15 kpc, something which is not seen in Figures |8] or |9l 
We therefore tentatively conclude that the upper boundary 
in Figure [121 is real. 

In Figure[l3]we have very crudely adjusted the K-hand 
bulge magnitude to the B-band by using B — JsT = 4, so as to 
place the bulges in the same diagram as elliptical galaxies. 
While they roughly overlap with the dwarf elliptical galax- 
ies, at a given magnitude they scatter to smaller sizes. This 
is somewhat reminiscent of their behaviour in the magni- 
tude - Sersic index diagram, where, for a given magnitude, 
the bulges scatter to smaller Sersic indices than the dwarf 
elliptical galaxies do (Graham 2001a, his Fig. 14). 



6 SUMMARY 

Most bright galaxies in the Universe today are spiral galax- 
ies. Loveday (1996) has estimated that two-thirds of galaxies 



brighter than AIv = — 18 mag are spiral galaxies and Lintott 
et al. (2008) report that the elliptical to spiral ratio is 0.57 
for galaxies with Mr < — 19 mag; the actual number of disc 
(i.e. spiral plus lenticular) galaxies is greater still. Discs are 
well known to contain dust which biases our measurements 
of their physical parameters and also that of their embed- 
ded bulges. Matching observational data from over 10,000 
galaxies (Allen et al. 2006) with radiative transfer measure- 
ments from sophisticated, three-dimensional models of disc 
galaxies with dust (Popescu et al. 2000; Tub's et al. 2004), 
Driver et al. (2007) has recently provided a highly useful 
calibration of these models. From this, equations to correct 
the observed scalelength and surface brightness for inclina- 
tion and dust are derived and presented here (equations [1] 
and [2}. In the optical bands, dust can bias the measured 
disc scalelength and surface brightness by 40 per cent and 
~0.5 mag arcsec"^, respectively, while the disc magnitude 
can be in error by up to 1 magnitude. Corrections for dust 
are therefore quite significant. We remark that programs to 
measure the evolution of disc galaxy parameters may need 
to additionally consider the evolution of their dust. 

These attenuation-inclination corrections enabled us to 
provide a modern, quantitative measure of the intrinsic 
physical properties of disc galaxies in the local {z < 0.03- 
0.04) Universe. We have done this by converting the ob- 
served properties of disc galaxies into their intrinsic (not 
simply face-on) dust free values. This is desirable because 
comparison between theory and observations, which is used 
to facilitate our understanding of galaxy formation and evo- 
lution, is hampered when one's observational measurements 
are compromised by the attenuating effects of dust. 

Using a sample of galaxies which have been carefully fit- 
ted with Sersic R^^" bulges plus exponential discs, this pa- 
per presents intrinsic structural properties for disc galaxies 
as a function of galaxy type. This work improves upon mea- 
surements obtained using R^^"^ models for every bulge, and 
provides inclination- and dust-corrected disc scalelengths, 
central surface brightnesses and magnitudes in optical and 
near-infrared passbands. 

In addition, we present inclination- and dust-corrected 
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bulge-to-disc flux ratios in several passbands, and K-hand 
bulge-to-disc size ratios, as a function of morphological type. 
For a given galaxy type, the logarithm of the _B-band bulge- 
to-disc flux ratios from Simien & de Vaucouleurs (1986) are 
as much as ~0.5 dex larger than the values obtained here. 
This is predominantly due to the best-fitting Sersic i?^/" 
bulge models having n < 4, while Simien & de Vaucouleurs 
force fit an i?^''* model to every bulge. In general, the me- 
dian dust-corrected bulge-to-total fiux (mass) ratios are less 
than 1/3 (1/2), which may provide a useful new restraint on 
some simulations which tend to produce larger ratios. Given 
the observed range of bulge Sersic indices from ^^0.5 to ^^4, 
we advocate the use of Prugniel & Simien's (1997) density 
model which was developed to match deprojected Sersic pro- 
files, rather than only matching a deprojected i?^''* profile. 

Rather than only dividing the galaxies into early- and 
late- type bins based upon their colour and/or concentration 
index, catalogued Hubble types were available for all of our 
galaxies. This has also enabled us to perform our analysis 
on a sample of disc galaxies, rather than an unknown blend 
of early- and late-type galaxies which have similar concen- 
tration indices at low magnitudes. However no attempt has 
been made to measure volume densities given the ad-hoc 
galaxy selection criteria. Although we note that in Driver et 
al. (2007, 2008), which contains expressions to correct bulge 
and disc magnitudes in a range of optical and near-infrared 
passbands, including the Sloan Digital Sky Survey (SDSS; 
York et al. 2000) filter set, we have already done this for the 
luminosity density of bulges and discs. 

We have presented updated (surface brightness)-size 
diagrams for bulges and discs and investigated the size- 
luminosity relations. We reveal that the discs of late-type 
disc galaxies occupy a broad region in these diagrams and 
are not well represented by linear relations. In contrast, the 
discs of early-type disc galaxies appear to reside over a nar- 
rower territory, defining the bright envelope to these distri- 
butions. We have additionally presented the size-luminosity 
relation for bulges and elliptical galaxies, emphasizing that 
a linear relation is inappropriate. The bulges of disc galax- 
ies are shown to overlap with, but additionally scatter to 
smaller sizes than, the distribution of similarly luminous el- 
liptical galaxies. 

To conclude, while the amount of dust in galaxy discs 
is roughly one thousand times less than the stellar mass, 
it can have a large influence on the observed properties of 
galaxies. In this study we have corrected for the influence 
of dust in order to obtain the intrinsic stellar properties of 
galaxies. Such characteristic physical properties, and indeed 
their bivariate distributions and possible scaling relations, 
are not only interesting in their own right but provide valu- 
able z = constraints for models of galaxy evolution. 
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